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Using the equation of motion technique for Green's functions we derive the exact solution of the 
boson fermion model in the atomic limit. Both (fermion and boson) subsystems are characterised 
by the effective three level excitation spectra. We compute the spectral weights of these states and 
analyse them in detail with respect to all possible parameters. 

Although in the atomic limit there is no true phase transition, we notice that upon decreasing 
temperature some pairing correlations start to appear. Their intensity is found to be proportional 
to the depleted amount of the fermion nonbonding state. We notice that pairing correlations behave 
in a fashion observed for the optimally doped and underdoped high T c superconductors. We try to 
identify which parameter of the boson fermion model can possibly correspond to the actual doping 
level. This study clarifies the origin of pairing correlations within the boson fermion model and may 
elucidate how to apply it for interpretation of experimental data. 



Boson fermion (BF) model describes a system com- 
posed of the narrow band electrons or holes (fcrmions) 
which coexist and interact with the local pairs (hard- 
core bosons) of, for example, bipolaronic origin [1]. The 
BF model has been recently intensively studied by var- 
ious methods, such as: the standard mean field theory 
[1] , the perturbative procedure with respect to the boson 
fermion coupling [2] , perturbative expansion with respect 
to the kinetic hopping [3] , the dynamical mean field pro- 
cedure [4], the continuous canonical transformation [5], 
etc. Apart of studying the mechanism responsible for 
superconductivity, there have been also investigated the 
many-body effects which, above T c , lead to an appear- 
ance of fermion pairs without their long range coherence. 
Indeed, three independent procedures [2,4,5] gave unam- 
biguous arguments for the precursor effects, out of which 
a pseudogap is the most transparent one. 

The pseudogap feature gradually builds up upon lower- 
ing temperature. It is observed in a temperature regime 
T* > T > T c , with both characteristic temperatures 
T* and T c depending on the BF model parameters. 
Absence of the long range coherence between pairs is 
caused by quantum fluctuations of the order parameter 
(cijCif) = Xi^ i ■ I n general, it is hard to distinguish be- 
tween the amplitude Xi and phase <pi fluctuations because 
they are convoluted. Intuitively one may expect that 
phase fluctuations would dominate for a dilute concen- 
tration of paired fcrmions, while in the opposite limit the 
amplitude fluctuations take over. Some analysis along 
this line was recently discussed in Ref. [6]. Fluctuation 
effects were also studied for the 2 dimensional (isotropic 
and anisotropic) BF model by Micnas et al [7] using the 
Kosterlitz Thouless theory. Authors reported a notice- 
able splitting between T* and T c which considerably in- 
creased for increasing population of the paired fermions. 
This result supports the above mentioned reasoning. 

In this brief report we show that already on a level of 
the zero-dimensional (atomic limit) physics there is some 



evidence for pairing correlations which gradually increase 
in strength upon lowering temperature. We study such 
effect on a basis of the rigorous solution of the BF model 
in the atomic limit. 

In our previous paper [3] we have investigated some as- 
pects of the atomic limit solution. The effective fermion 
spectrum was determined there by a direct diagonalisa- 
tion of the Hilbert space. In a current work we rederive 
the exact solution using the equation of motion technique 
[8] for Green's functions. Advantage of this method is 
that it gives the spectral weights for the eigenstates ex- 
pressed in terms of the corresponding correlation func- 
tions. Of course, diagonalisation and Green's function 
method are equivalent and complementary to each other. 

Hamiltonian of the boson fermion model can be writ- 
ten as H = X^j, ct ^j c 1,ct c .j> + where stands 
for the hopping integral and the local part Hi is given by 

[1] 

H = e ^2 c t<r c i,o- + E b\bi 

a 

+ 9 ( b i c lAd + 6 Kf c u) ■ (!) 

We use here standard notations for the second quantisa- 
tion operators of fermion c,* )CT , c\ a and hard core boson 

bi, b\ fields. Site energies are correspondingly expressed 
as Eq — Sf — /U and Eg = As — 2/i where a common chem- 
ical potential \i ensures conservation of the total charge 
concentration n to t = (%b\bi + J2 a c \ <r c i,<^j- Fermion and 
boson fields are coupled through the exchange interaction 
gbicl »ct . + h.c. which can transform a fermion pair into 
a hard core boson and vice versa. 

In the strict atomic limit = one needs a solution 
of only the local part (1). Let us notice that the hard 
core boson operators obey, in general, the spin \ alge- 
bra, characterised by the following commutation rules 
[bub]] = 8 l0 {l - 2b\b t ) and [b h b s ] = = [&+,&}]. For 
the same site i = j (which is relevant in the atomic 
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limit) they simply reduce to the anticommutation rela- —iQ(t)( Ai(t),A\(0) ) = J duie iuJt ((Ai; Aj))^. 
tions [9]. We can thus construct the fermionic Green's 

function ({Ai;A\)) u both for fcrmions A, = c tlJ and According to the equation of motion [8] w((A; B)) u = 

for hard-core bosons A { = b u where we introduced ({A B}) + (([A, H] ; B)) u wc find the following set of cou- 

the Fourier transform of the retarded Green's function pl°d equations 

(w-eo)((ci lT ;ct T )) w = l + g{{bic\ A ;c\ A )) w , (2) 
(w + eo-^)<<6ict i ;4 t )) u ,= fl <<(< i -nf) 2 c i , T ;4 T )) u , (3) 
(w - E ) «(nft - nf ) 2 c iiT ; 4 T )) W - ((n& - nff) + <7«&4 a ; 4 T » W , (4) 

where n F a = c] CT Cj )(T and nf = 6-6^. After some algebraic calculations we determine that these three functions read 



((ci, T ;4 T )) w = ! + 



W - £0 



(w - e ) (w + s - E ) - g 2 



It is convenient to rewrite the single particle Green's function in the following way 

Z F 



uj + e Q - E 



«Ci,T;ct T )) w 
Z F 



UJ -£ 



+ (1 - Z') 



+ 

UJ — S + UJ — S- 



1 - <« x - nff) 
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V = 1 - U 
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(5) 
(6) 

(7) 

(8) 
(9) 

(10) 
(11) 



Another set of coupled equations to determine the hard core boson propagator ((bf, b\)) u involves the following Green's 
functions 



(w-Eo) = l + g{{ci d a,T,b\)) u , 

{uj - 2s ) ((q,i c *,t; & !»^ = 2 (c ia ci, T 6}) +g{{b 1 ;b^) u - 5^((ct <T c i , ff 6 i ; 6t)) w , 

<7 

{uj - E ) £« c t 6j)) w = « T + n^) . 



(12) 
(13) 

(14) 



In analogy to (8) we present the explicit form of the single particle Green's function as 



+ (1-Z i 



uj — E 

ryB II F F \2\ 

Z = (K,t ~ n i,l) > ' 
E± = s± + s . 



uj — E + uj — E- 



(15) 

(16) 
(17) 



The single particle propagators (8) and (15) are both 
characterised by a three pole structure. One of the poles 
is a remnant of the free nonbonding state (eo for fermions 
and Eo for hard core bosons). The other two poles 
(s± and E±) correspond to the bonding and antibond- 



ing states which arise due to the boson fermion interac- 
tion. Hamiltonian (1) is no longer diagonal in the oc- 
cupation representation \n F ,n F ;n B ) because two eigen- 
vectors contain admixture of | T>l>0) and |0, 0; 1) [3]. 
Loosely speaking, an ability of the system to fluctuate 
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between these two states is a measure of pairing correla- 
tions (we mean the correlations in time, because in the 
atomic limit there exist no spatial correlations). 

Let us inspect in some detail the spectral weight Zp 
of the nonbonding fermions' state. From (9) we see 
that Z F is depleted from unity by {{nf ^ — nf ) 2 ) = 

( n fl) ( n F) ~~ (^ n i \ n f) '■ -^ means that propagation (in 
time) of the free fermion (with spin a =\) occurs unless: 
(a) there exists another fermion on the same site with 
the opposite spin and simultaneously no hard-core bo- 
son is present there, (b) there is boson while J, fermion is 
absent. Disappearance of the nonbonding state depends 
thus on fermion and boson concentrations. Role of other 
factors, such as for example temperature, is less evident 
at this point. 

Spectral weight of hard core boson nonbonding state 
is given by 



Y) = {nt)-2{nr r ) , 



(18) 



where nf 



fermions on site i, while n^ atr counts only the doubly 
occupied fermion states nf" r = c\ »ct .Cj^c^f. The hard 
core boson can safely exist in a free (nonbonding) state 
when there are only single fermions present on the same 
site. The more fermions are paired, the less spectral 
weight is left for a free hard core boson. 

We can express the spectral weights Z F and Z B 
explicitly via the concentrations n F = J2a( c l <j c i,<?) > 
n B = (b\bi) and through such parameters as tempera- 
ture T and Ab- From a general relation [8] (AB) = 
— i J du)f(u))Im&g((B; A)) u+iri we obtain 



Z* = 



n F - p/( £+ )+ M 2 /( £ _)] 
f(e )-[v\f(e + )+^f(e-)} ' 

n B ~ [u*f(E + ) + v 2 f(E-)] 
f(E )-[ U 2f(E + )+v*f(E_)} ' 



(19) 
(20) 



where f(x) = [e xf3 + l] 1 is the Fermi Dirac distri- 
bution and /3 = l/k B T. These quantities can be 
computed also from the diagonalized Hamiltonian us- 
ing the Lehmann representation. They are found to 



be [3] Z F = [1 
(6 



1 + 2e 



/© 
+ 



_ -0eo _|_ e -/3(eo+-E ) _|_ e -/3(2e +-Eo) j 
-Peo + 2e~ l3( - eo+E ° S> + e -P( 2 zo+E ) 

e -P E + + e~P E - is the partition function) and Z B = 
[2e- f3£o + 2e-^( £o+i30 )] /Q. These expressions are of 
course identical with (19,20). 

We explored numerically variation of the spectral 
weights Z F , Z B versus temperature T and Ab for sev- 
eral fixed charge concentrations n to t = n F + 2n B . From 
our analysis it turns out that the most sensitive In- 
dependence of these quantities occurs for £ + Eq = 
when n tot — 2. One can show that 

2 "" (21) 



J \n to 



3 + 



sh/3- v /(A B /2)2+^ 
cosh(/3A B /6) 



- 7 B 



which at high temperature approach the asymptotic 
value limT^ooZ F ' B _„ = 0.5, while for T — ► diminish 
to zero. Figure 1 illustrates this behaviour. 
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FIG. 1. Spectral weight of the nonbonding state of the 
fermion and hard core boson subsystems for total charge con- 
centration ntot = 2. Main suppresion of the spectral weight of 
the nonbonding state occurs near T* (pointed by the arrows) 
and depends on the parameter As. 



In any other case the spectral weights Z F , Z B may 
not vanish in the ground state. They vary within a nar- 
rower regime signalling that interaction effect is then less 
efficient as compared to the case n tot = 2. Figure 2 
shows the spectral weights Z F ' B as functions of n tot for 
A^/2 = Sf. For fermions we notice that away of n tot — 2 
the spectral weight Z F increases and becomes less de- 
pendent on temperature. In the extreme dilute region 
Z F — ► 1. As far as Z B is concerned it follows the be- 
haviour of Z F only in a close vicinity of n to t — 2. Going 
away from such case the nonbonding spectral weight Z B 
decreases as a direct consequence of the relations (16,18). 
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FIG. 2. (a) Spectral weight of the nonbonding°fermion 
state as a function of total charge concentration ntot per site. 
All finite temperature values are situated within the shaded 
area in the figure, (b) Spectral weight of the nonbonding 
state in the hard core boson subsystem for several tempera- 
tures as indicated. Both figures (a) and (b) are obtained for 
As/2 = s f . 

Parameter Ag has rather a minor effect on both spec- 
tral weights, it mainly affects their temperature variation 
similarly to what is shown in figure 1. In order to charac- 
terise the temperature dependence of Z F we define char- 
acteristic temperature T*(ntot, As) [3] which is an in- 
flexion point d 2 Z F {T*)/dT 2 = 0. Roughly speaking, the 
spectral weight Z F starts to decrease when temperature 
drops below T* . From higher dimensional studies of the 
BF model in the symmetric case (n to t — 2, = 2sf) 
[3,4] it is known that suppression of the nonbonding state 
below T* is accompanied by an appearance of the pseu- 
dogap structure. Apart of the symmetric case there is 
not enough evidence that such relation remains valid. 

Reduction of the nonbonding state spectral weight Z F 
for temperatures near and below T* is closely related to 
appearance of the pairing- type correlations. To prove 
this let us consider the Green's function ((6^c| ^ ; ct 
given in equation (6) which yields the following correla- 
tion function 

(clcl lbi )= 9 (l-Z F ) fi£ f + Z f £ { _ £ - ] ■ (22) 

Let us recall that on a level of the mean field theory [1,7] 
the superconducting order parameter is given as 

k 

where £k = \J (ik — A 1 ) 2 + |5(^)I 2 an d £k denotes a dis- 
persion of itinerant fermions. In the atomic limit the 
order parameters are on average equal zero {c\^c\^j = 
= (hi). We can think of a finite value (22) as a result 
of fluctuating pairing correlations. Magnitude of pair- 
ing correlations vanishes at high temperatures while, for 



temperatures T < T* , achieves the finite value propor- 
tional to the spectral weight depleted from the nonbond- 
ing state (1 — Z F ). Figure 3 illustrates the temperature 
dependence of pairing correlations for several values of 
As at the fixed charge concentration n to t = 2. Mag- 
nitude of (cj |cj turns out to be proportional to the 
mean field value of T C MF which proves their close relation. 
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FIG. 3. The pairing correlation function (c? jC? ^6,) in- 
duced in the atomic limit by the boson fermion coupling g. 
The upper inset shows a mean field value T C MF for a fixed 
ratio g/D — 0.1, where D denotes the fermion bandwidth. 
The lower inset is the ground state value of the correlation 
function. The main figure and the insets were obtained for 
ntot = 2. 



The boson fermion model is claimed by some authors 
[1,6,7,10] to capture key aspects of the theory for high 
temperature superconductors (HTSC). In realistic de- 
scription of the HTSC materials one must however con- 
sider their anisotropic dim = 2 + 5 structure. Pairing 
correlations discussed here for the atomic limit would in 
higher dimensions lead to: (a) formation of fermion pairs 
at T p , and (b) at T c < T p to their long range coher- 
ence, establishing the superconductivity (with T c ^ 0). 
What remains to be studied for the realistic dim = 2 + 5 
systems is a pseudogap region of the incoherent fermion 
pairs T c < T < T p . We hope that the exact solution of 
the BF model discussed here for the atomic limit may 
help in such future investigations. 
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